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Recently, in tlie context of f{R) modified theories of gravity, a new type of model has been 
proposed where one directly couples the scalar curvature to matter. As any model in f{R) theory, 
there are certain conditions which have to be satisfied in order to ensure that the model is viable and 
physically meaningful. In this paper, one considers this new class of models with curvature-matter 
coupling and study them from the point of view of the energy conditions and of their stability under 
the Dolgov-Kawasaki criterion. 
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I. INTRODUCTION 

Recently, there has been a revival of interest in grav- 
ity models where the linear scalar curvature term of the 
Einstein-Hilbert action is replaced by a function, f{R), of 
the scalar curvature. This interest is due to the possibil- 
ity of accounting for the late time accelerated expansion 
of the universe without the need of explicit additional de- 
grees of freedom such as dark energy and also to possibly 
replace dark matter at galactic level and beyond (see, for 
instance, Refs. [l|, 0| for reviews) . 

In this paper, one examines models where, besides the 
extension in the geometric part of the action, there is also 
a direct coupling of a new function of the Ricci scalar to 
the matter Lagrangian density 3 (see also Ref. ^ for a 
recent review): 



(1) 



where fc is a coupling constant, fi{R) (with 1 = 1,2) are 
arbitrary functions of the curvature scalar R and £,„ is 
the Lagrangian density corresponding to matter. One 
sets k = 1 hereafter. The second function, f2{R), is 
usually considered to have the following form. 



/2(i?) = 1 + A^2(i?) 



(2) 



where A is a constant and ip2 is another function of R. 
Notice that this non-minimal curvature-matter coupling 
has been previously considered in order to address the 
issue of the accelerated expansion of the universe Q . 
Varying action |T]) with respect to the metric g^^j/ yields 
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the field equations 



(/l + 2£m/2)i?piy — -^figfiv — ^fififl + 2£m/2) — /2'7)iiy , 

(3) 

where the prime denotes the differentiation with respect 
to R, A^^ = V^V^ - g^^D and □ .g'^'^V^V^. The 
matter energy-momentum tensor is defined as usual 



T — 



2 5{^/^C„i) 



(4) 
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This theory exhibits interesting features. The most 
intriguing is the fact that the motion of a point-like test 
particle is non-geodesic due to the appearance of an extra 
force in the theory Actually, this issue has been a 
lively point of discussion in the literature [1, 0, B S ■ 
These modified gravity models have also been examined 
from the point of view of their impact on stellar stability 

m- 

In this work one analyzes these new modified models 
of gravity with non-minimal coupling between curvature 
and matter from the point of view of the energy condi- 
tions and their stability under the Dolgov-Kawasaki cri- 
terion. 

As one shall see, the energy conditions allow one to 
establish under which conditions gravity remains attrac- 
tive, along with the demands that the energy density is 
positive and cannot flow faster that light, in the context 
of the new modified models, while the Dolgov-Kawasaki 
stability criterion helps one ruling out some classes of 
modified models p^ . 

This paper is organized as follows. In section [TT] one 
investigates the energy conditions suitable for the mod- 
ified gravity models with non-minimal curvature-matter 
coupling. One deduces an equivalent energy-momentum 
tensor, along with an energy density, pressure and the 
gravitational coupling and thus obtain the energy con- 
ditions for the modified gravity models. At the end of 
section ini one illustrates the energy conditions by apply- 
ing them to a fairly broad class of models. In section Hill 
one studies the generalized Dolgov-Kawasaki criterion for 
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some modified gravity models and examine their viabil- 
ity in what respects their stability. Section IIVI contains 
our conclusions. 



II. ENERGY CONDITIONS 

A. The Raychaudhuri Equation 

The origin of the null energy condition (NEC) and of 
the strong energy conditions (SEC) is the Raychaudhuri 
equation together with the requirement that the gravity 
is attractive for a space-time manifold endowed with a 
metric g^i/ (see e.g. Ref. for a discussion). From 
the Raychaudhuri equation one has that the temporal 
variation of the expansion of a congruence defined by the 
vector field is given for the case of a congruence of 
timelike geodesies by 



d9 
d7 



In the case of a congruence of null geodesies defined by 
the vector field k^, the Raychaudhuri equation is given 
by 



de 
d7 



1 



a^,a^''' + u^.tu'^'' - R^^k'^k'' , (6) 



where Rf^u, 0, a^i, and lo^i, are the Ricci tensor, the ex- 
pansion parameter, the shear and the rotation associated 
to the congruence, respectively. 

It is important to realize that the Raychaudhuri equa- 
tion is purely geometric and independent of the grav- 
ity theory under consideration. The connection with the 
gravity theory comes from the fact that, in order to re- 
late the expansion variation with the energy-momentum 
tensor, one needs the field equations to obtain the Ricci 
tensor. Thus, through the combination of the field equa- 
tions and the Raychaudhuri equation, one can set physi- 
cal conditions for the energy-momentum tensor. The re- 
quirement that gravity is attractive imposes constraints 
on the energy-momentum tensors and establishes which 
ones are compatible. Of course, this requirement may not 
hold at all instances. Indeed, a repulsive interaction is 
what is needed to avoid singularities as well as to achieve 
inflationary conditions, and to account the observed ac- 
celerated expansion of the universe. 

Since cr^iyCT^'' > 0, one has, from Eqs. ([5]) and ©, that 
the conditions for gravity to remain attractive {dO/dr < 
0) are given by 



Rf.^u^'u" > 
R^^kf'k" > 



SEC 
NEC 



(7) 
(8) 



for any hypersurface of orthogonal congruences (w^^ — 
0). 

For the sake of completeness, one can easily see that 
if one inserts Einstein's field equations into Eq. ([7]) one 



gets, 



T, 



T 



(9) 



where T^"^ is the energy-momentum tensor, T is its trace 
and one has used units where 8ttG = c = 1 (which shall 
be adopted hereafter). If one considers a perfect fluid 
with energy density p and pressure p, 



(10) 



then Eq. ^ turns into the well known SEC for Einstein's 
theory. 



p + 3p>0. 



(11) 



In the same fashion, if one inserts the Einstein's field 
equations into Eq. ([S]), one obtains 



'^ni^k^ k ^ , 



(12) 



which, after assuming Eq. (jlOp . turns into the familiar 
NEC of general relativity 



p + p>0. 



B. Effective Energy-Momentum Tensor 



(13) 



In what follows, one generalizes to the gravity model 
with non-minimal curvature-matter coupling the proce- 
dure developed in Ref. [l3| for f{R) theories. 

Rewriting Eq. ([3]) in order to make explicit the Ein- 
stein tensor, G^t^, one gets 



Gfj,i, — 



f2 



f'l ^" 2/2-^71 



(14) 



where the effective energy-momentum tensor T^i, has 
been defined as 

f,. = m - iill^R^ 5..+-^ A,. (/{ + 2C^f^) 

(15) 
ig 

(16) 



^From Eq. (|14p one can define an effective coupling 

A + 2i2m/2 

The field equations can be written as 

G^. = k (f^. + Tp.) . (17) 

Thus, in order to keep gravity attractive, besides con- 
ditions Eqs. (I?])-®, k has to be positive from which 
follows the additional condition 



/2 



> 0. 



(18) 
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It is important to realize that this condition is indepen- 
dent from the ones derived fr'om the Raychaudhuri equa- 
tion (Eqs. dll) and ([8])). While the latter is derived di- 
rectly from geometric principles, the former is related 
with the fairly natural definition of an effective gravita- 
tional coupling. 

The effective energy-momentum tensor defined by Eq. 
(If 5p can be written in the form of a perfect fluid (Eq. 
(fTO| ) if one defines an effective energy density and an 
effective pressure. However, given the presence of the 
higher order derivatives in Eq. p^ . in order to proceed 
it is necessary to specify the metric of the space-time 
manifold of interest. Since one is interested in cosmolog- 
ical applications, the Robert son- Walker (RW) metric is a 
natural choice. In what follows, one considers the homo- 
geneous and isotropic flat RW metric with the signature 



dt^ ~ a^{t)dsl 



(19) 



where ds^ contains the spacial part of the metric and a{t) 
is the scale factor. 

Using this metric the higher order derivative term is 
given by 

A^,/i(ii', Cra) = (V^ V, - g^,,U)h{R, Cm) 

= {d^d^ - gf,udodo)h - (r°^, + g^^3H)doh 

(20) 

where h{R,£m) is a generic function of R and H = 
d/a is the Hubble expansion parameter and F^^ = ad6^jy 
(with /1, 1/ ^ 0) are the components of the affine connec- 
tion. For the flat RW metric one has R — -~6{H^ + f)- 
Thus, the effective energy density is given by 

- _ 1 / /l fj + 2£m/2 p \ Q tj /f + 2^™/2 u /r,-|\ 

^ - 2 U h " h ^''^ 



while the effective pressure is given by 

^ I fi f j + 2/3m/2 

/2 



R 



■ f" + 9r f" f" -\- 9r f" ■ 

{R + 2HR) + ^^"-^^ + ^1 +^^'»/2 ^ (22) 
where the dot refers to derivative with respect to time. 

C. Strong and Null Energy Conditions 

^From the above results, it is straightforward to obtain 
the energy conditions. From Eq. the Ricci tensor 

can be written as 

= +T^. -^^^-(T-fT)^ . (23) 

Inserting this into Eq. ([7]) one obtains for SEC 

A(T^. + T^.)M^M"-^fc(T + T) >0. (24) 



Using Eq. (fT8|) and assuming a perfect fluid as a matter 
source, one gets 



/2 



+3{R + HR) 



h 

/i + 2/3,„/2 
J2 



R 



fill ,cyr rill 

h 

In the same way, inserting Eq. (|23|) into Eq. 
obtains the NEC 



(25) 



(26) 



and hence 



p + p+{R- HR) 



^•'-mJ2 d2 



h 



/2 



i?^ > 0. 
(27) 

It is important to realize that these results are ob- 
tained directly from the Raychaudhuri equation. How- 
ever, equivalent results could have been obtained from 
the transformations p ^ p + p and p p + p of the SEC 
and NEC of general relativity. 

Naturally, setting /2 = 1 in Eqs. ^ and ^ one 
recovers the results of Ref. [l3|. Setting fi — R and 
/2 = 1 one recovers SEC and NEC of general relativity 
(Eqs. ([ni) and (HSl)). 

It is interesting to note that if one imposes that gravity 
is repulsive, the signs of the inequalities ([24|) and ((26|) 
would be the same, since the signs of the conditions ([7]), 
(HI) and (flS)) would change, and thus the above results 
for the SEC and the NEC would remain unaltered. 



D. Dominant and Weak Energy Conditions 

For the derivation of the DEC and WEC, one can ex- 
tend the results obtained above and apply the transfor- 
mations p ^ p + p and p p + p to the DEC and WEC 
of general relativity since it is natural to assume that the 
energy conditions will not be violated when one changes 
from the Jordan frame to the Einstein frame (see, for 
instance, Ref. [11]). 

Thus, one gets for the DEC 



p-p + 



./2 



fi + 2£„t/2 

/2 



R 



{R + 5HR) 



2-Cm/2 



fill 
Tl 



f2 



-|^^i?2>0 (28) 



and for the WEC 

1 ^h_f{ + 2Cmn 
J 2 



P 



R] - 3H 



fl + 2£m/2' 



R>0. 



2 V/2 ./2 J f2 

(29) 

As one may easily check, for /2 = 1 one obtains the 
results of Ref. and if one further sets /i = R, Eqs. 
and (P5|) yield p — p>0 and p > 0, as expected. 
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E. Energy Conditions for a Class of Models 

To get some insight on the meaning of the above energy 
conditions, one apphes then to a specific type of models 
where /i 2 are given by 



fi{R) = R + eR^, 



(30) 



where 'P2{R) is defined by Eq. ([2]). Assuming a i 
RW space-time, Eq. ([T9| . the energy conditions can 
written in the following way. 



be 



2XC 

n 



l + A|i?|^ 



> b (31) 



where a, b and ck„^„i depend on the energy condition un- 
der study and one defines e = (— l)"e and A = (— 1)'"A 
due to the fact that for a RW metric one has R < 0. For 
the SEC, one finds 



,SEC 



= -1, b 



SEC 



= -{p + 3p), 



a. 



SEC 



-n[l + 3(n - l)(ij + HR)R-^ + 
S{n^l){n-2)R-^R^]. 

For the NEC, one obtains 

^NEC _ n i^NEC 



(32) 



(33) 



0, 



3(n- 1)HR-^R 



(35) 



a^^^ = -n{n - 1) {R - HR)R-^ + {n ~ 2)R-^R^ 
For the DEC, one has 

a^^^^l, b^^^^-ip-p), (36) 

a^E'^ = n{l + {n- l)(i? + 5HR)R-^ + 

{n-l){n-2)R-^if). (37) 

Finally, for the WEC, one gets 



(38) 



(39) 



Writing the derivatives of R in terms of the decelera- 
tion (q), jerk (j) and snap (s) parameters. 



1 a 

H'' a 



1 'a 1 "a , , 

H-^ a a 



(34) one obtains for a„. 



-1 + — [n - 1) - — TTi^(n - l){n - 2) 



2{q-lY 



2{q-lf 



(41) 



NEC 



n{n — 1) 



6(g-l)2 



6(g-l)3 



{n-2) 



(42) 



g^ + 3g-H5j + 5-4 ^ , {j-q-2)\ . 
1 TTT^ (n - 1) H 77 TT^l"- ^ ^ 2) 



6iq-iy 



6(g- 1)3 



(43) 



2"2(T3^("-^) 
I 



(44) 



Given these definitions, the study of all the energy 
conditions can be performed by satisfying the inequal- 
ity ([3T|) . In Table U it is presented the energy conditions 
for models of type Eq. ([50]) . 

One realizes that all the energy conditions depend on 
the coupling constants of the model and on the space- 
time under consideration. The balance between the cou- 



plings and the space-time parameters and their evolution 
along the history of the universe will determine whether 
the energy conditions are satisfied or not. Setting the 
estimated values of g, j and s from observation, one can 
impose bounds on the coupling constants of the gravity 
model (see e.g. Refs. [13, ESEa] for discussions). 

When applying for models like Eq. (|5D|) . the condition 
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e > 


e < 


A > 


2XCtii a — GiLrL^-^^\ri — m 


2A£^ ^ aipin — m 
le — a„i 1 1 




1\\\Cjn ckti— alpln — m 
e — OLjri ' ' 


2 AZ^Tn a — Qj^l^n — m 
e — 0:7,1 1 ' 


1 - \x\\Rr < 


2|A|£^rji ^- q:ti— aipin — m 





TABLE I: SEC, NEC, DEC and WEC expressed by Eq. ([SlJ with 6 = and where it is assumed that > 0. For < 
one just needs to change the direction of the inequahties. 



for attractive gravity (AG), Eq. (fTS]) . gives, once again, 
an inequality like Eq. (pij) with 



= 1 : 



.AG 



6H^{l-q) 



(45) 



(46) 



Thus, the results depicted in Table U also stand for the 
condition that ensures that gravity remains attractive. 



III. GENERALIZED DOLGOV-KAWASAKI 
INSTABILITY 

Modified gravity theories may be subjected to instabil- 
ities. These are due to the fact that the field equations 
([3]) are of fourth order and their trace give a dynamical 
equation for R. Such instabilities do not occur in general 
relativity since there the field equations are of second or- 
der and their trace give an algebraic equation for R. In 
order for the dynamical field R to be stable its "mass" 
must be positive, a requirement usually referred to as 
Dolgov-Kawasaki stability criterion. 

The Dolgov-Kawasaki criterion was studied in the con- 
text of f{R) theories in Ref. [17] and generalized for the 
models under study here, i.e. models whose action is 
given by Eq. ([T]), in Ref. [l^l and can be expressed as 



/{'(i?) + 2A/:„.^^'(i?) >0 



(47) 



where <f2{R) is defined by Eq. ^ and fi{R) is written 
as fi{R) = R + eipi{R), with e being a constant. Notice 
that Eq. ([TT]) is not quite the one presented in Ref. 
given that there the Lagrangian density. Cm, was not 
considered. 

For models such as the ones defined by Eq. (|5D|). the 
inequality (|T7|) becomes 

en{n - l)|i?|"-2 + 2XC„rm{m - l)|i?r" > , (48) 

where e and A are defined by, 

. [(-l)"e, ifi?<0 ^ [(-1)™A, ifi?<0 
^ ^ I e, if i? > ' ^1 A, if i? > ' 



Notice that this inequality can also be written as Eq. 
(ED) for. 



0, b 



DK 



0, 



(50) 



a. 



DK 



n{n — 1) 
36^4(1 -g): 



r(l + A|i?r 



(51) 



Thus, the results presented in Table HI for the cases 
with l-hAli?!™ > 0, stand for Eq. ^ too. However, it is 
important to realize that although the energy conditions 
and the Dolgov-Kawasaki criterion yield the same type of 
inequalities, they are independent from each other. From 
TableU one can see that, if the factors are the same, 
the various conditions (the four energy conditions, the 
condition for attractive gravity and the Dolgov-Kawasaki 
criterion) are degenerate. 

For most of the cases, the viability of the model with 
respect to the Dolgov-Kawasaki instability criterion will 
depend not only on the value of the constants e and A 
(which may be further constrained from the solar sys- 
tem observations [loj|). but also on the space-time metric 
under consideration. 

^From the Dolgov-Kawasaki criterion one has 



n{n — 1) 
m(m — 1) ' 



(52) 



and thus the stability criterion for models defined by Eq. 
(|30)) do not depend on the values of j or s, and the space- 
time dependence is present only through R. Further- 
more, since a^^ /a^^ > 0, one can see that models for 
which A > and e > 0, the Dolgov-Kawasaki criterion is 
always satisfied independently of the values of e, A and R 
(note that as < 0, the direction of the inequalities 

in Table [T] are reversed). 

For models where m = n the inequality Eq. (I48p im- 
plies, e -|- 2XCjn > 0. 

Another type of models that is interesting to consider 
are the ones where /i,2(-R) are given by 



/i(i?) = ^a„i?" 

n=l 
k' 

m— 1 



(53) 



For these models, the inequality Eq. ([47]) becomes 



(49) J2 " l)i?""^ + 2A£™ J2 bmm{m - l)R"'-^ > 



n=2 



m=2 



(54) 

The stability conditions for models with k,k' = 2,3 
are presented in Table |TT1 
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k' = 2 


fe' = 3 




aa + 2\Cmb2 > 


as > 


^ ^ a2+2A£^fa2 
— 3a3 


a -i < 


— 3|a3| 


k = 3 


XCmbi > 


^ — 6AC„63 


as + 2\Crabz > 


p ^ a2+2A£T7^62 
- 3{a3+2A£„i,3) 


XLmh < 


- 6|A£mb3l 


as + 2\Crabz < 


p ^ a2+2A£TTi&2 
- 3|a3+2A£„i,3| 



TABLE IL Stability conditions for models defined by Eq. (|54p . where it is assumed that 7? > 0. 



Once again, one sees that the stability conditions for 
most of the models are metric dependent. This fact is 
in strong contrast with what happens for /(i?) theo- 
ries where the stability can be ensured independently of 
the space-time. This difference can be relevant and it 
brings about further constraints on /(i?) theories with 
a matter-curvature coupling since it restricts the Ricci 
scalar, which might turn out to be unphysical. 

IV. CONCLUSIONS AND OUTLOOK 

In this work, the well known energy conditions and 
the Dolgov-Kawasaki stability criterion, for the modified 
theories of gravity described by action ([T]) were examined. 
In section |ll] it was shown that the generalized SEC and 
NEC can be obtained directly from the SEC and NEC of 
general relativity through the transformations p — > p + p 
and "p — *■ p + p, where the quantities p and p contain 
the physics of the higher-order terms in R. The same 



procedure can be applied to the DEC and the WEC. As 
expected, all generalized energy conditions for a suitable 
energy-momentum tensor are strongly dependent on the 
geometry. Conditions to keep the effective gravitational 
coupling positive and gravity attractive (Eq. (|18p ) were 
also obtained. The encountered energy conditions have 
been applied to a specific kind of models and shown that 
all the energy conditions together with the positiveness 
o f the effective gravitational constant and the Dolgov- 
Kawasaki criterion can be expressed by the same type of 
inequalities. 

In section IIIII the Dolgov-Kawasaki stability criterion 
was used to test some modified gravity models. One finds 
that the stability of a specific model will depend not only 
on its couplings but also on the space-time geometry in 
question as depicted in Table [TTl 

We stress that our method and in particular Eqs. (|18p . 
(PS)) . ((77)) . (^5]) and ^ are quite general and can be used 
to study the physical implications of any f{R) model with 
non-minimal curvature-matter coupling 
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